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Self-similarity of clusters of galaxies and the Lx — T relation 
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Abstract. In this paper based on ROSAT/PSPC data we investigate the emission measure profiles of a sample 
of hot clusters of galaxies (kT > 3.5 keV) in order to explain the differences between observed and theoretically 
predicted Lx~T relation. Looking at the form of the emission measure profiles as well as their normalizations 
we find clear indication that indeed the profiles have similar shapes once scaled to the virial radius, however, 
the normalization of the profiles shows a strong temperature dependence. We introduce a M gas -T relation with 
the dependence M gas oc T 1 ' 94 . This relationship explains the observed Ly^-T relation and reduces the scatter in 
the scaled profiles by a factor of 2 when compared to the classical scaling. We interpret this finding as strong 
indication that the Af gas -T relation in clusters deviates from classical scaling. 
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1. Introduction 

The simplest, purely gravitational, models of structure 
formation predict that galaxy clusters are self-similar in 
shape and that scaling laws relate each physical pro perty 
to the cluster total mass M and re dshift z (Kaiser 
Navarro et al 



1997; Teyssier et al. 1997; Eke et al 



1986 



1998 



Bryan & Norman 1998 ). Self similarity applies to both the 
dark matter component and the hot intra-cluster medium 
(ICM). The virial theorem then yields the well known scal- 
ing relations: 



M s oc T 3/2 (l + z)- 3/2 
Rs oc 7^/2(1+^-3/2. 



(1) 

(2) 



where Ms is the total mass in the sphere of radius Rs 
corresponding to the overall over density S and T is the 
cluster temperature. Assuming a constant gas mass frac- 
tion, the total gas mass then scales as M gas oc T 3 / 2 and 
the X-ray luminosity as T 2 . 

The comparison between observed clusters and theo- 
retical predictions gives us insight into the physics that 
governs their formation and evolution. Vikhlinin et al. 
1999| recently found R oc T°- 57±0 04 , which is close to the 



predicted scaling law. Both the gas density and tempera- 
ture profiles of hot clusters (T > 4keV) do show regularity 
(Markevitch et al. |1998| Neumann & Arnaud [1999| , here- 
after paper I; Vikhlinin et al. 1999[ ). The shapes of vari- 
ous clusters, once the radius is scaled to the virial radius, 



look remarkably similar outside the cooling flow region, 
supporting the existence of an universal underlying dark 
matter profile. 

However, clusters also deviate from the predicted scal- 
ing laws. The most remarkable deviation is the observed 
Lx-T relation. Different authors, studying different clus- 
ter samples, found Lx oc T 2 ' 9 (Arnaud & Evrard 1999; 
Markevitch |1998j ; Allen & Fabian [1998D with a low dis- 
persion, ruling out the theoretical relation Lx oc T 2 . This 
difference between the slopes of the observed and theo- 
retical relations can be explained if i) the cluster shape 
depends on the ICM temperature and/or ii) the gas mass 
temperature relation deviates from predictions (Arnaud & 
Evrard |l999| ; Ponman et al. |l999| ). Non gravitational ef- 
fects, like pre-heating by early galactic winds, could affect 
the gas distribution shape (e.g Loewenstein 200C; Tozzi 
& Norman 2001) and have been proposed to explain the 
steepening of the Lx~T relation (Kaiser 1991; Evrard & 
Henry 1991). However the observed structural similarity 
of hot clusters suggest that pre-heating plays an impor- 
tant role only for cool (T < 3.5 keV) clusters (see also 
Ponman et al. h999L Lloyd-Davis et al. 12000). 
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In order to investigate further the observed Lx~T re- 
lation, we examine in this paper the emission measure 
profiles of a sample of nearby hot clusters, for which the 
effect of pre-heating is a priori minimal. The shape and 
normalization of the emission measure profiles are sen- 
sitive to both the cluster internal structure and the to- 
tal gas mass. Comparing the profiles of clusters at differ- 
ent temperatures can provide more detailed information 
about the cause of the discrepancy between the predicted 
and observed Lx~T relations. In this paper we adopt 
Hq = 50km/sec/Mpc and £l m = 1. 



2 



D.M. Neumann, M. Arnaud: Self-similarity of clusters of galaxies and the Lx — T relation 



2. The data 

The emission measure profiles are from our previous work 
(Paper I), where the cluster sample and imagery data pro- 
cessing are described in detail. The cluster sampl e com- 
prises Abell clusters (Abell, Corwin & Olowin 1989 ) in the 
redshift range 0.04 < z < 0.06, which were observed in 
pointing mode with the ROSAT PSPC and provide good 
statistics. Here we only consider the clusters in the sub- 
sample for which accurate temperature measurements ex- 
ist from the literature (with kT > 3.7keV). The emission 
measure profiles are deduced from the observed surface 
brightness profiles via Eq.3 in paper I. 

For the Lx~T relation, we adopt here and in the fol- 
lowing the relation derived by Arnaud 



L x otT 



2.88 



Evrard ( ^99g ): 
(3) 



which is in good agre ement with other works (Markevitch 
19981 ; Allen & Fabian |l99S| ). 



3. Combining self-similarity and the observed 
Lx — T-relation 

3.1. Dependence of X-ray luminosity and emission 
measure profile on gas mass and ICM structure 

Let us consider a cluster of extent R and temperature T. 
Its X-ray luminosity depends on its total gas mass (within 
radius R) and its internal structure as: 



r ^as 



A(T)Q(T) 



(4) 



(see also Arnaud & Evrard 1999 ) where A(T) is the cooling 
function, which we will assume to be T 1 / 2 . Q(T) is a form 
factor which only depends on the shape of the gas den- 
sity distribution (p gas (r)). Q{T) is equal to (p 2 gas ) / (p g as) 2 , 
where the brackets denote the a verage over the whole clus- 
ter (see also Arnaud & Evrard 1999 ). If clusters are self- 
similar in shape, Q(T) should be constant. 

The emission measure along the line-of-sight at radius 
r is defined as: 



EM (r) oc 



n 2 {x)xdx 



Vx 2 



M 2 

gas 



F(r/R,T) 



(5) 



F(r/R,T) is again a dimensionless form function, which 
should be independent of temperature if clusters are self- 
similar. According to the theoretical scaling laws, EM 
should scale as T 1 / 2 . 



3.2. Does the ICM shape depend on temperature? 

In Paper I, we already quantified the structural variation 
in the cluster sample. We derived the logarithmic ICM 
density gradient a by fitting a /3-model to the data, 



a(r/R) = - 



_ d\ogn g (r/R) q J3 

d\og(r/R) ~ ''l+(r c /ry 2 ■ 



(G) 



and showed that the dispersion of a is less than 20% at any 
scaled radius. The a values, for two different scaled radii, 
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Fig. 1. The a parameter (logarithmic density gradient) 
versus temperature. The crosses show a at a scaled radius 
of r/R=0.3, the circles at r/R=0.6. 

are plotted versus cluster temperature in Fig.|l|. There is no 
apparent correlation between a and T , i.e no systematic 
variation of shape with temperature. 

Fig.|| shows the scaled emission measure profiles of the 
clusters, EM(r/R) oc EM(r) /T 1 / 2 (upper panel), and 
their corresponding relative dispersion (lower panel). The 
scaling radius was fixed to R = T2oo with the normaliza- 
tion from Evrard et al. (1996). Due to detection limits 



we only consider radii r/R < 0.5, our sample becoming 
incomplete at larger radii. The dispersion of the profiles 
is significant (40% standard deviation) but stays constant 
with radius. This again indicates that the profiles are es- 
sentially parallel (self similarity of form) but points to- 
wards a possible problem with the normalization scaling 
factor. These features are also evident to the eye in the 
upper panel. 

We will thus assume in the following a common shape 
for all the clusters in the sample (i.e. Q{T) and F(r/R, T) 
are constant) and study possible deviation from the theo- 
retical M gas -T and R-T relations. 

3.3. Temperature dependence of the gas mass 

We now consider: 



R oc T x and 



oc 



(7) 



From Eq. ||, Eq. [| and Eq. |?], the x and y slopes must 
be related by 3x — 2y — 0.62 to account for the observed 
slope of the Lx — T relation. 

If we assume that x — 0.5 , corresponding to the clas- 
sical R-T relationQ we derive y = 0.44 or Af gas oc T 194 . 
From Eq. |E|, the emission measure scales in this case as 
EM oc T 1 - 38 , instead of EM oc T 1 / 2 . The scaled emis- 
sion measure profiles, using this new scaling, are shown 
in Fig.||. Note that these rescaled profiles are simply the 

1 Another formal solution is y = 0, which implies x = 0.21. 
This is in contr adicti on with the R — T relation found by 
Vikhlinin et al. ( 1999 ). Therefore we discard this possibility 
in the following. 
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Fig. 2. Upper panel: Scaled emission measure profiles 
using the theoretical scaling laws (arbitrary units). For 
A3562 we used kT = 5.2keV a new measurement obtained 
by Ettori et al. 2000. Lower panel: The corresponding 



standard deviation divided by the mean value as a func- 
tion of scaled radius. 



scaled profiles, EM, displayed FigJ|, divided by T 88 . 
Remarkably, this re-normalization does reduce the disper- 
sion of the profiles, by a factor of 2, when compared to the 
classical scaling. 

To study more directly the deviatiqnjrom the theo- 
retical temperature scaling, we display EM , estimated at 
r/R = 0.3, versus T in Fig.j§ To extend the upper limit of 
the sample temperature range, we also considered A2163 
(kT ~ 13keV), one of the hottest clusters observed so far 



(Elbaz et al. 1995 ). A clear correlation between EM and 
T is observed. The relationship required to explain the 
observed Lx — T relation (EM oc T ' 88 , full line) matches 
the data points. The dashed line in Fig, fi represents a 
crude best fit to the data (Press et al. 1993) ignoring er- 
ror bars. This fit gives EM oc T°", reasonably close to 
the required relationship. The error bars on EM and T are 
highly correlated in a non trivial way and we did not try 
to estimate error bars on the correlation. At fixed scaled 
radius, EM is derived from the EM value, measured at a 
physical radius which depends on T via the R-T relation. 
This can be an important source of error, not included 
here, since EM varies rapidly with radius. The relatively 
large dispersion at low temperatures (Fig.|J) is possibly 
due to the larger number of clusters in this domain and 
our underestimate of the true errors. 



4. Discussion 



Recently Mohr et al. ( 199E ) estimated the gas mass within 
^500 for a sample of 45 galaxy clusters. Their derived 

1 ' s±018 , is steeper than the 



Fig. 3. Upper panel: Scaled emission measure pro- 
files (arbitrary units) assuming M gas oc 7 11 - 94 ( x = 0.5, 
y = 0.44). For A3562 we used the kT measurement by 
Ettori et al. |2000[ The dynamic range is the same as 
in Fig j^. Lower panel: The full line shows the corre- 
sponding standard deviation divided by the mean value 
for all clusters. The dotted line shows the relative disper- 
sion without A780, which hosts a strong cooling flow and 
which increases the dispersion considerably at small radii. 




Fig. 4. The scaled emission measure of the clusters (clas- 
sical scaling law) versus temperature at r/R = 0.3 (about 
1 Mpc). The error bars for the temperature are 90% c.l. 
and for the emission measure la. The full line represents 
EM oc T - 88 and the dashed line EM oc T°" (see also 
text). 



theoretical expectation and in excellent agreement with 
our results. Our results are also consistent with the rela- 
tion, M gas OC 



M gas -T relation 



M gas cx T 1 



r i.7i±0.i3 ) d erive d b y Vikhlinin et al. ( |l999| ) 
after fitting a power law to the X-ray emission of the outer 
parts of clusters. Although this general agreement was not 
unexpected, since all these studies are based on ROSAT 
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surface brightness data, it is not entirely trivial. Our evi- 
dence of a steeper than predicted M gas -T relation relies on 
structural and scaling properties of the emission measure 
profiles, guided by the independently measured Lx~T re- 
lation, while the analysis in previous works is more global, 
based on the integrated density profiles. 

Mohr & Evrard (1997) found a tight correlation be- 
tween the X-ray isophotal radius and the temperature 
with Ri oc T - 93 ^ - 11 . This relation is steeper than the 
theoretical expectation, R\ oc T 2 / 3 , for self similar clusters 
with = 2/3 (Mohr et al. |199S|). They interpret their find- 



ing as evidence for structural variation with temperature. 
However, our study indicates that this steepening (for 
clusters with kT > 4 keV) originates from systematic devi- 
ation of the profile normalization from the theoretical scal- 
ing with T. For a /3-model , the X-ray surface brightness 
profiles at large radii varies as I(r) oc EMo(r/R c ) 1 ~ ei3 , 
where R c is the core radius and EMq the central emission 
measure. For self similar clusters R c oc R$ oc T 1 / 2 and 
Ri scales as EM^T 1 / 2 (for = 2/3). If EM scales as 
T 1 - 38 , instead of T 1 / 2 , as derived here, we ex pect t hat Ri 
scales as T 096 , as found by Mohr & Evrard Qt997D. 



The steepening of the M gas -T relation has to be ex- 
plained physically. It can arise from i) a variation of the 
gas mass fraction / gas with temperature and/or ii) a de- 
viation of the total mass versus T relation from the theo- 
retical scaling law. 

If the modeling of dark matter clustering is correct 
and the ICM is in hydrostatic equilibrium, the latter pos- 
sibility would imply that the X-ray mean temperature 
is a biased indicator of the virial temperature. For in- 
stance relativistic particles creating radio halos could con- 
tribute to the ICM pressure. If the strengths or quantity of 
these particles i s a fu nction of temperature as indicated by 
Colafrancesco ( 1999|) it might explain the deviation from 
the M — T relation. Departure from isothermality could 
also bias the temperature Tx measured from single tem- 
perature fits from global cluster X-ray spectra. Numerical 
simulations (Mathiesen & Evrard 2001 Bialek et al. 2001 ) 
suggest a slight steepening of the M — Tx relation with 
M oc T x e ~ 17 . New X-ray observatories such as XMM or 
Chandra will be able to investigate this possibility in de- 
tail. 

The effect of pre-heating on the M gas -T relation was 



studied recently by Bialek et al. ( 2001 ). Numerical simula- 
tions with initial entropy level consistent with the observa- 
tions of Lloyd-Davis et al. (2000), predict a steepening of 



the M gas -T relation, typically M gas oc T 1,90±0 ' 07 , which 
is coherent with our results. The M-T relation obtained 
from simulations with and without pre-heating is similar. 
However, in pre-heating models, the M gas -T relation be- 
comes steeper than the M — T relation, due to systematic 
variations of / gas with T. However in the work by Bialek 



et al. ( |200l|) only 7 clusters have temperatures above 3.5 



keV, with a large scatter in the / gas ~T relation (figure 
13, their work). Larger samples of simulated hot clusters 
and more detailed studies of their internal shape are nec- 
essary to see if pre-heating can contribute to the bias in 



the M gas -T relation via / gas variation, without introduc- 
ing structural variations, which is inconsistent with the 
observations. 



5. Conclusion 

Studying the emission measure profiles of a sample of clus- 
ters (kT > 3.5keV) we found a strong indication that the 
Mgas-T relation of hot clusters deviates from theoretical 
expectations. Our results suggest M gas oc T 1 - 94 . This re- 
lationship together with the fact that clusters show self- 
similar profiles can explain the observed Lx~T relation, 
which also differs from theoretical expectations. 

This result is in agreement with previous work such 



as Mohr et al. (1999) or Vikhlinin et al. (1999), based 



on the direct measurement of M gas in clusters. However, 
our study is more sensitive to deviations from the pre- 
dicted M gas -T relation since the emission measure is, with 
EM oc Mg as , strongly coupled to the gas mass content. 
Observations from XMM-Newton or Chandra of the tem- 
perature structure are required to assess whether the total 
mass also diverges from the predicted scalings. 
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